Long-Range Effects in Layered Spin Structures 
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We study theoretically layered spin systems where long-range dipolar interactions play a relevant 
role. By choosing a specific sample shape, we are able to reduce the complex Hamiltonian of the 
system to that of a much simpler coupled rotator model with short-range and mean-field interactions. 
This latter model has been studied in the past because of its interesting dynamical and statistical 
properties related to exotic features of long-range interactions. It is suggested that experiments could 
be conducted such that within a specific temperature range the presence of long-range interactions 
crucially affect the behavior of the system. 
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Many of the forces that we see in the universe have a 
long-range nature where in d dimensions a pairwise inter- 
action potential decays as V(r) ~ i/ r d + s with — d < s < 
0. Examples include gravitational interactions, Coulomb 
and magnetic forces. Nowadays there is a vast litera- 
ture which describes various exotic nonlinear and sta- 
tistical properties of many-body systems with long-range 
interactions such as disagreement of predictions of micro- 
canonical and canonical ensembles, negative specific heat 
and temperature jumps, etc. However, most of these 
studies are purely theoretical and there have been only 
few suggestions on how to test experimentally all these 
peculiarities. One example is the astrophysical observa- 
tion of negative specific heat [2] which is an outcome of 
the truly long range nature of the gravitational force. In 
addition, it has been suggested that systems composed 
of a small number of particles can show negative specific 
heat Q due to non-additivity even when the interaction 
is short range. This has been verified experimentally in 
nuclear collisions [3, atomic sodium clusters Q and in 
molecular clusters [6|. However no experimental test of 
these predictions has been carried out for a laboratory 
system with long range interactions. 

This Letter aims at proposing testable effects of dipo- 
lar magnetic forces, whose long-range nature is a conse- 
quence of the cubic decay law of the potential V(r) ~ 
1/r 3 . This law results in a strong dependence of the 
dipolar energy on the sample shape (see e.g. Ref. 0]). 
However, in ordinary magnetic systems, dipolar energy 
is about a thousand times smaller than Heisenberg ex- 
change interactions between nearby spins. Therefore, in 
most cases the role of long-range forces is to introduce 
some anisotropy which determines the ordering direction 
in the magnetic sample. On the other hand, in nuclear 
magnets, where magnetic order is fully defined by dipolar 
interaction, one has to go down to nano Kelvin temper- 
atures to observe ordering 





FIG. 1: Left graph: schematic arrangement of spin s — 1/2 
Cu +2 ions in a (C v H2 V +iN H3)2CuCln layered compound. 
Right graph: Suggested form of the sample which allows the 
observation of the spontaneous flips of the magnetization vec- 
tor (thick arrows at the middle) which are predicted by Hamil- 
tonian ([7]) to which the microscopic model ([l]) effectively re- 
duces. 



In this Letter, we propose to examine more 
closely magnetic layered spin structures (e.g. 
{C v H2 V +iNH 3 ) 2 CuCh [|) in which the effective 
magnetic interaction between electronic spins is predom- 
inantly dipolar. In particular, we examine rod shaped 
layered spin structures, whose microscopic Hamiltonian 
can be effectively reduced to that of a one dimensional 
coupled rotator model with both nearest neighbor 
coupling and a dominant mean-field interaction term 
resulting from the dipolar forces. We suggest that these 
compounds provide a system in which exotic phenomena 
that characterize coupled rotator models with both 
short and long-range couplings [10, EH could readily be 
observed. Here, we propose to verify experimentally the 
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presence of all-to-all mean-field couplings by monitoring 
the time-scales of spontaneous magnetization flips below 
the magnetic transition. Similar effects of magnetization 
reversals could in principle ta ke p lace in systems with 
only short-range interactions [12j. However, here we 
propose to explore the experimental conditions in 
which collective reversals are driven by the presence of 
long-range interactions, and therefore their existence 
becomes strongly dependent on the shape of the sample. 
Moreover, the average reversal time is expected to 
grow as the exponential of the volume of the sample as 
opposed to the case of short range interactions where 
it grows only as the exponential of the surface area 
[12j. In addition, magnetization reversals due to long 
range interactions are expected to be sharper, since they 
involve at once the entire sample. 

The magnetic arrangement of the class of compounds 
(C u H2v+iN H 3 ) 2 CuCl4 is schematically given in Fig. 1 
(see for more details Refs. [HI, Hi]]). The system is com- 
posed of ferromagnetic layers with strong intralayer in- 
teractions, W, and a weak coupling, w, between the lay- 
ers which is either ferromagnetic for v = 1 or antifer- 
romagnetic for v > 1. This allows us to consider the 
three-dimensional spin system at temperatures or ener- 



gies well below W as a quasi-one dimensional ferromag- 
netic or anti-ferromagnetic spin chain consisting of clas- 
sical spins, with each spin in the chain representing a 
whole layer. The forces in the system are provided by 
the exchange interactions between the effective spins in 
the chain (short-range forces) and the dipolar interac- 
tions among all spins (long-range forces). The interlayer 
exchange interaction turns out to be comparable with 
the dipolar interactions making these systems excellent 
candidates for considering long-range effects, even in the 
case of a small number of layers. Note that, as expected, 
the shape of the sample is rather important in systems 
with dipolar interactions. A particularly interesting case 
is that of ferromagnetic interlayer coupling (y = 1) with 
a sample shape for which dipolar forces favor ferromag- 
netic ordering such as in a rod cut along the layer planes 
(rod along axis z in Fig. 1). This case will be discussed 
in detail below. 

These systems have been modelled previously includ- 
ing only the anisotropic contribution of the dipolar forces 
(see e.g. Ref. [l3j]) while the long-range character of these 
forces was neglected. At this stage we include fully the 
dipolar interactions and the Hamiltonian reads 



I,<i,j> I,<i,j> H^Jj 



where the first sum represents the intralayer exchange in- 
teraction (indices i and j refer to nearest neighbor spins 
within the same layer, while the indices I and J number 
the layers), [Ib is Bohr's magneton and r is the vector 
between the sites of the spins sjj and sjj. The param- 
eters £ and 77 (< 1) yield biaxial hard axes anisotropics 
along the x and y axes, respectively. They are a result of 
the crystallographic forces and take into account out-of 
plane and in-plane anisotropies (note that the main con- 
tribution to the in-plane anisotropy comes from dipolar 
forces themselves, see below). The second sum stands for 
the interlayer exchange interactions (with i and j now re- 
ferring to nearest neighbor spins in adjacent layers) and 
the last term describes dipolar interactions among all 
spins. According to [UHl W ^ 10 4 w in the compounds 
{C\H^N H^CuCl^. At low temperatures all spins in a 
single layer are ordered ferromagnetically and therefore 
the spin vector sjj could be considered as independent 



of the index i, sj = sjj, and this represents the spin of 
a whole layer. This is justified only if one works below 
the ordering temperature of a single layer, which approx- 
imately coincides with the intralayer exchange constant 
W. Under the additional condition nw < W (n is the 
number of spins in a single layer), in a certain tempera- 
ture range each layer will be ordered ferromagnetically, 
while the transition to 3D ordering will be strongly af- 
fected by the existence of long-range dipolar forces (which 
are comparable with the short range interlayer exchange 
w). We consider the thermodynamic properties of the 
system under such conditions. Applying the ordinary 
procedure to calculate the dipolar sum, one can divide it 
into a short-range contribution (restricting ourselves to 
consider the interaction between nearest neighbors) and 
a long-range one 0]. Then ([TJ) can be rewritten in the 
following one-dimensional representation (see Ref. [la ]) 
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where B x = 4W(1 - £) and B y = 4W(1 - ??) define the 
hard axis anisotropics along a; and y, respectively (the 
out-of plane anisotropy being much larger than the in- 
plane one B x 3> By), N is the number of layers, tu ex is 
an effective exchange constant between the layers con- 
sisting of the sum of the exchange constant w between 
the neighboring spins of the different layers and the con- 
tribution of dipolar forces between the same spins. Thus 
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Where q a and r a are the distances between nearest neigh- 
boring spins within the same layers and different layers, 
respectively. In the case of the considered compound 
q a = 5.25-lCr 8 cmandr a = 9.97-lCT 8 cm [H]. Finally the 



last two mean-field terms come from the long-range part 
of dipolar forces. Moreover, u>m = (iir/3)(2(i%/vo) and 
vq is the volume of the unit cell of the lattice. Note that 
the nearest neighbor part of dipolar forces favors antifcr- 
romagnetic ordering of the spins in neighboring layers, 
while the long-range component of the forces is ferromag- 
netic along z. The measured value for the out-of plane 
anisotropy is B x — 240mK, while the in-plane anisotropy 
constant B y is less than 5mK. The exchange constants 
is oj ex = 4.8mK (see e.g. Ref. 
Thus, we have: 



13j), and u>m = 16mK. 



B x /uj m = 15, By I Lo M rj 0.2, 



5 /% - 0.3. (4) 



As is evident, the out-of plane anisotropic term is associ- 
ated with a much larger energy scale than all the others 
and therefore we have neglected in |(3J) all the terms which 
include the s x component of the spins which emerge from 
dipolar or exchange interlayer forces. 



We now consider the torque equation 



dsj 



sj x Hj 



(5) 



where Hj — —dTL/dsj is an effective magnetic field act- 
ing on the spin sj. Noting that \sj\ = 1/2, it is natural to 
adopt the definitions Sj = 2s j, 
and 51 
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y/l-(S5) 2 sm0j 



< j — yl — {Sj) 2 cos 9j which yield the equations of 
motion in terms of angular variable 9j and of Sj. Ob- 
serving that B x is the largest energy scale in the Hamil- 
tonian 
as follows 



one can simplify the equations of motion for Sj 



~~dT 



nB x Sj 



(6) 



and approximate -Jl — (Sj) 2 ~ 1 ■ Then, substituting Sj 
from the last equation into Hamiltonian @, we obtain 
the following effective Hamiltonian, which includes only 
the angles 0j's, and their time derivatives 
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where the dimensionless time is introduced via the trans- 
formation t — > tny/B x ujM- The above model describes a 
system of coupled rotators with nearest neighbor and all- 
to-all mean-field interactions. 

In analogy with the similar system considered in 
Ref. (loj , we can derive the thermodynamic and dynami- 
cal properties of the model ([7]). It follows that, in the en- 
ergy range —0.3 < e < e c = 0.376, two maximal entropy 

N 

states with opposite magnetization m = ( cos 9j)/N 

J=i 

are present and, in principle, for finite N, the system 
can flip between these states. Above the critical en- 
ergy e c the system is characterized by a single maxi- 
mum entropy state at m = 0. At e c a second order 
phase transition of the mean-field type is present. The 
lower energy limit derives from the fact that, in the ther- 
modynamic limit (JV — > oo), the energy e satisfies the 
condition e > —m 2 /2 — aj ex /u>M- It follows that below 



e = -Wa/wjf — —0.3 there is an ergodicity breaking 
region [11], where the system retains its magnetization 
direction for infinitely long time. If the system energy 
is slightly above the ergodicity breaking limit, magneti- 
zation flips are very rare. Simulating the system using 
torque equation ([5]) does not reveal any flip for observable 
time scales, i.e. the system behaves like in the ergodicity 
breaking region. 

Examining the dynamics of the one-dimensional spin 
chain ^ at e = 0.2, sufficiently above the ergodicity 
breaking region, the time evolution of the magnetization 
displays the expected behavior (see upper graph in Fig. 
2). We then compute numerically the probability distri- 
bution function (PDF) of m, P(m), and present the data 
of (1/N) log(P(m)) as circles in the inset. These are com- 
pared (modulo a vertical normalization shift) with the 
analytical curve for the entropy derived for Hamiltonian 
([7]). The agreement is impressively good, considering the 
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FIG. 2: Main plots: Numerical simulations of the one di- 
mensional spin chain model @, (JSJ) with different energies: 
e — 0.2 (upper graph) and e = 0.4 (lower graph), i.e. below 
and above the second order phase transition energy e c = 0.376 
predicted for Hamiltonian (J7J). The number of layers in the 
simulations is taken N = 100. Insets show the corresponding 
entropy curves (solid lines) and data obtained from the PDF 
of the magnetization (open circles). 



absent. 

In summary, we have studied theoretically the layered 
spin structure (C v H2v+iN H^)2CuCli and showed that, 
under certain conditions, it could be modeled as a sys- 
tem of coupled rotators with both short and mean-field 
interactions. In experiments, this could be confirmed by 
the presence of magnetization flips for v = 1 (ferromag- 
netic interlayer short-range interactions). When v > 1 
the short-range interlayer coupling is antiferromagnetic, 
while each layer is still ordered ferromagnetically. For 
an appropriate choice of the short versus long-range cou- 
pling (which could be precisely controlled by changing 
the length of rod shaped sample), the latter compounds 
could serve as a laboratory system for which various ex- 
otic phenomena, such as ensemble inequivalence, nega- 
tive specific heat, temperature jumps, etc. 2, 1(| , char- 
acterizing systems with long-range interactions, could be 
observed. 
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numerous approximations we have made when deriving 
Hamiltonian (J7J) from ©. 

If now one further increases the system energy above 
the critical one, e c , magnetization fluctuates around zero 
and the corresponding PDF is presented in the lower 
graph of Fig. 2, together with the analytical prediction, 
again in very good agreement. 

Let us conclude with a few remarks that are relevant 
for the experimental implementation. Let us recall that, 
throughout this paper, we have considered a rod shaped 
sample. Moreover, as stated above, the reduction of the 
Hamiltonian to a one-dimensional rotator model holds 
only if nw < W. Thus, for the compound under consid- 
eration, n should be less than 10000. Hence, in a real 
experiment, one should examine samples of typical size 
n = 20 x 500 within a single layer and, e.g., N — 20 lay- 
ers in total. In contrast, in spherical samples with about 
the same number of spins per layer, long-range effects 
will be negligible and only short-range forces are left in 
the effective one-dimensional model: no phase transition 
will be present and, therefore, magnetization flips will be 
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